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ABSTRACT 


The aim of this letter is to discuss the virtual identity of two recent tidal theories: our creep tide theory and one Maxwell 
model recently developed. It includes the discussion of the basic equations of the theories, which, in both cases, include 
an elastic and an anelastic component, and shows that the basic equations of the two theories are equivalent and differ 
by only a numerical factor in the anelastic tide. It also includes a discussion of the lags: the lag of the full tide (geodetic), 
dominated by the elastic component, and the phase of the anelastic tide. In rotating rocky bodies not trapped in a 
spin-orbit resonance {e.g., the Earth) the geodetic lag is close to zero and the phase of the semidiurnal argument in the 
anelastic tide is close to 90 degrees. The results obtained from combining tidal solutions from satellite tracking data and 
from the Topex/Poseidon satellite altimeter data are extended to determine the phase of the semi-diurnal argument in 
the Earth’s anelastic tide as cro = 89.80 ± 0.05 degrees. 


1. Introduction 


the external body; i?2,o is the eccentricity function 


In his celebrated work on the secular change of elements 
due to tides, Darwin (1880) introduced ad hoc lags in the 
potential of a tidally deformed body. These lags were small 
and proportional to the frequencies of the delayed terms of 
the potential. Over the years, the fact that these two at¬ 
tributes were hypotheses introduced by Darwin has been 
forgotten. During the past century, some models with con¬ 
stant lags were attempted, but a serious discussion of the 
proportionality to frequency hypothesis had to wait for the 
work of Efroimsky and Lainey (2007) showing that in the 
case of terrestrial planets and planetary satellites, an in¬ 
verse power law would be a better choice. But the lags 
continued to be small and the smallness of the tidal lags 
remained as a principle carved in stone. 

The contradiction between being small and following an 
inverse power law dominated discussions in the past decade 
and even served to argue the impossibility of the inverse 
power law. New results in contradiction with that postulate 
came from the creep tide theory (Ferraz-Mello, 2012, 2013). 
One of the intermediary results of the creep tide theory is 
the shape of the body deformation due to the anelastic tide, 
which, after the transient phase (i.e. for jt ^ 1), is dom¬ 
inated by the semidiurnal component. On the equatorial 
plane, the semidiurnal anelastic tide is 


Can = -j^Re^pE2fl COS (Tq COs{2ip - 2£ - (Tq) 


( 1 ) 


where Zp is the prolateness of the ideal Jeans spheroid rep¬ 
resenting the static tide due to an external body M at the 
distance a (semi-major axis); ip is the longitude of one point 
on the surface of the body; Re is the mean equatorial ra¬ 
dius; e is the orbital eccentricity; i is the mean longitude of 
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V 
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( 2 ) 


where 7 is the relaxation factor (also known as critical fre¬ 
quency) and V = 2D — 2n is the semidiurnal frequency (D is 
the rotation velocity of the body and n is the orbital mean- 
motion). We remember that in the creep tide theory, ag is 
a fully determined constant introduced by the integration 
of the creep equation and not an ad hoc plugged lag. 

The maximum of Can is reached when 2p — 2i — ag = 0, 
i.e. the angle between the vertex of the point where the 
height is maximum to the sub-M point is (To/ 2. In the case 
of giant planets, the critical frequency 7 is in the range 
10 — 100 s“^ (see Ferraz-Mello, 2013). The frequencies in¬ 
volved in the tide (rotation, mean-motion) are in the range 
10“® — 10“"^ s“^. Therefore, in this case v y and so 
(To ^ 0 and the anelastic tide highest point remains almost 
aligned with the mean direction of the tide raising body M. 
However, in the case of rocky planets and planetary satel¬ 
lites, 7 is in the range 10“® — 10“^ s“^. Then y v and 
CTo will approach 90 degrees. Besides, in intermediary cases 
approaching synchronization, ctq may take any values in the 
interval (0,7r/2) and may lead to significantly large geodetic 
lags. 

This last result contradicts the assumption of small¬ 
ness of the lag postulated by Darwin, and some current 
beliefs. In the discussion of the creep tide theory, Ferraz- 
Mello (2013) stated that “this result is in contradiction with 
the observations. For instance, the observed geodetic lag of 
the Earth’s body semidiurnal tide is very small (0.16 ±0.09 
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degrees cf. Ray et al. 1996)". This comment refers to the 
comparison of the geodetic lag with the lag of the anelastic 
creep tide and is not valid out of that context. In reality, 
as discussed in the next section, it is not in contradiction 
with the observations, but just with one particular inter¬ 
pretation of them. The value of o-q close to 90 degrees is 
consistent with the observations and cannot be considered 
as a setback of the creep tide theory as stated by Correia 
et al. (2014). 


2. The observed value of the Earth's tidal lag 

We start this section quoting a statement found in Zschau 
(1978): “Measurements of tidal gravity variations at the 
Earth’s surface, as well as precise observations of the tidal 
effect on satellite orbits have not yet revealed reliable re¬ 
sults on imperfectly elastic body tides on Earth". In 2001, 
referring to earlier attempts to determine the energy dis¬ 
sipation in the Earth’s body, Ray et al. (2001) said, “Un¬ 
fortunately, none of these early attempts to deduce k 2 from 
satellite tracking data was successful". Only by the end of 
the century, combining tidal solutions from satellite track¬ 
ing data and from Topex/Poseidon satellite altimeter data, 
Ray et al. (1996) managed to separate ocean and Earth 
tide signals and to determine the Earth body’s dissipation. 
Their more recent result (Ray et al. 2001) corresponds to a 
geodetic lag of the Earth’s body semidiurnal tide equal to 
£o = 0.20 ± 0.05 degrees. 

We recapitulate briefly how the lag was determined in 
the Ray et al. papers. They compare the Earth potential 
sensed by satellite tracking, which includes the contribu¬ 
tions of the Earth’s body, oceans and atmosphere, to the 
ocean component estimated from the Topex/Poseidon al¬ 
timeter data. The difference is due to the other components 
which are accordingly modeled. The anelastic part of the 
tidal potential of the Earth’s body is modeled using a clas¬ 
sical Darwinian model. Using the notations of Ferraz-Mello 
et al. (2008), the potential component corresponding to the 
lunar semidiurnal tide, in the planar approximation, is 


The angle (Tq is not an ad hoc lag, but a well-determined 
parameter that, in the case of the Earth, is very close to 90 
degrees. This equation can then be expanded to become 


SU 2 P = — 


2kfGmE? 
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0 E 2 .P cos (To sin( 2 (/j — 2i) sin ag ( 6 ) 


where we have neglected one term factored by cos^ ag ~ 
^ 0. We have the same equations as before, but 
instead of sineoj we have cosnosinno = |sin( 7 r — 2cto). 
When this is used in the Ray et al. equations, the result 
is (To = 89.80 ± 0.05 degreesQ The different approach for 
the introduction of the elastic tide in the creep tide theory 
does not affect the result because it was equally introduced 
in both tracking and altimeter data and only the difference 
between these two values actually matters in the determi¬ 
nation of the lag. 


3. The equations of the creep tide theory 

The creep tide theory of Ferraz-Mello (2013) is founded on 
one equation: the Newtonian creep differential equation 


C = -7(C-P) (7) 

where C.{g},0) is the height of the anelastic tid 43 at one 
point at the surface of the body, p{ip, 6) is the corresponding 
height of the ellipsoidal figure of equilibrium due to the 
joint action of tide and rotation, and 7 is a relaxation factor 
inversely proportional to the equivalent uniform viscosity of 
the body. The solution is the anelastic tide, which is then 
added to the elastic tide \p (where A is a free parameter 
related to the height of the tide) to give the final result. 

We can merge these two tide components into only one 
equation, li Z = ( -\- Xp, it is easy to see that 

Z-\-jZ = {I-\-X)jp-\-Xp ( 8 ) 

or. 


U 2 .g = - epEz.o cos( 2 y 7 - 2i - eg) (3) 

or, since eg 1 , we obtain the anelastic contribution: 


Z = Ce -I- Ap -I- 76 / pe^^dt. 


(9) 


In order to see that Eq. m is the equation of a Maxwell 
body, it is enough to replace p by the stress 


SU 2 .g = epE' 2,0 sin(2(p - 2i) sineo (4) 

where kf is the Love number for a homogeneous body and 
G the gravitational constant. The static term (correspond¬ 
ing to £0 = 0 ) does not need to be considered as it was 
already taken into account in the two potential data being 
compared. The anelastic part of the potential introduces 
sin£o in the Ray et al. equations, which, after elimination 
of the other effects (ocean loading and atmospheric tide and 
loading), is determined as £0 = 0.20 ± 0.05 degrees. 

However, we can proceed in a different way and use 
the creep tide theory instead of the Darwinian model. In 
that case, following Ferraz-Mello (2013), the component of 
the potential corresponding to the semidiurnal tide, in the 
planar approximation, is 


2kfGmR^ _ /n \ 

0U2,g = -- epE2,g cos (To COS(2(p -2£-(Tg). 
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(5) 


r = p-c = (i + A)p-z, ( 10 ) 

in the creep equation. Equation (I5|) then becomes 

z = (i + A)7r+Ar, (11) 

which is the constitutive equation of a Maxwell body (see 
Verhas, 1997). 

^ The comparison of the two results, valid when ag ~ 7r/2, gives 
Eo = 7r/2 — (jQ. 

^ In Ferraz-Mello (2013) ^ was the distance to the center of 
the body and the height of the tides were represented by S(^ 
and h (see Eqs. (55) and (58) in that paper). The change of 
the origin from the center of the body to one fixed reference 
level (e.g. a sphere with same volume as the body) simplifies 
notations when anelastic and elastic tides are composed and has 
no consequence in the study of the creep because both sides of 
Eq. ((TJ are invariant to such translation to the extent that p is 
also referred to the same level as C,. 
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3.1. Comparison to the equations of the Maxwell model 

In order to compare the equations of the creep tide theory to 
those of the Maxwell model used by Correia et al. (2014), we 
need to rewrite Eq. (9) of Correia et al. using as a variable 
the deformation C instead of the potential Vp. Proceeding 
exactly as was done in that paper, but keeping ^ in the 
right-hand side of their Eq. (6), we obtain 

C 4 C = “Ce H-Ce- (12) 

T T T 

If we substitute the two free parameters t and Te by those 
used in the creep tide theory using the equivalence formulas 

7=-, A=^, (13) 

T T 

we obtain, for the basic equation of Correia et al. (2014), 

C + 7C = 7Ce + Ad (14) 

whose solution is 

( = -h ACe + 7(1 - A)e-T'‘ J pe'^*dt. (15) 


Before comparing this result to that of the creep tide model, 
we have to discuss the meaning of C in both theories. In 
Correia et al. (2014), ( is the radial deformation of the free 
surface. In Eerraz-Mello (2013), C is also the deformation 
of the free surface, but before the inclusion of the elastic 
part in the theory. So, the ( of the Maxwell model is to be 
compared not to C of the creep equation, but to Z. We also 
note that = p. 

The comparison of Eqs. m and m shows that the 
two theories are virtually identical. The only difference be¬ 
tween the two theories is the numerical factor (1 — A), which 
appears in Eq. (1151) multiplying the part of the solution cor¬ 
responding to the anelastic tide. (We note that 0 < A < 1). 
This explains why the results of Correia et al. (2014) for 
the dissipation and the rotation look identical to those of 
Ferraz-Mello (2013, 2015). 

As an example, the empirical formulas relating the creep 
theory parameters to the dissipation parameter h^jQ of the 
classical theories are 

^•2- _ 1 7 ^ 

Q - 


in Ferraz-Mello et al. (2013) and 


Q 


fcy(l-A) 


71Z 


7^ -I- 


in Correia et al. (2014). The equivalence between these two 
equations was already shown by Eq. (90) in Correia et al. 
(2014). 


3.2. The Newtonian creep of the Maxwell model 

The inverse of the transformation used in Section |3] can 
be applied to Eq. m to give the equation of the creep 
embedded in the Maxwell model of Correia et al. If we 
denote by Can the radial deformation of the surface due to 
the anelastic component, that transformation is 


C — Can T ACe 


(16) 


and the result is 

Cln = -7(Can-(l-A)Ce). (17) 

This is a Newtonian creep which differs from that consid¬ 
ered in Eq. © only by the fact that here the stress is taken 
proportional to the distance to a different equilibrium sur¬ 
face, defined by (1 — A)p, instead of the surface of a Jeans 
ellipsoid, p, as in Eq. ©. 

The solution of Eq. dm is trivial. It gives, for the semid¬ 
iurnal anelastic tide of the Maxwell model, 

Can = ^-Reep£’2,o(l “ A) COS (To cos(2v? - 2£ - (To) (18) 

where (Tq is the same angle defined by Eq. 0. This means 
that the phase of the anelastic tide in the Maxwell model 
of Correia et al. (2014) is the same as that of the creep tide 
theory of Ferraz-Mello (2013). 


4. Conclusions and summary 

1. The dynamical tide can be decomposed into two parts: 
one elastic part, which corresponds to the perfect deforma¬ 
tion of the body under the tidal stress, and the so-called 
anelastic part, which corresponds to having the body per¬ 
manently adjusting its shape to follow the changing tidal 
potential. In the two theories discussed in this letter, Ferraz- 
Mello(2012, 2013) and Correia et al. (2014), these two 
parts are virtually the same. The only difference between 
them is that the solution given by Correia et al. for the 
anelastic tide is the same solution found in the creep tide 
theory of Ferraz-Mello multiplied by the numerical factor 
(1-Te/r) = (1 - A)_. 

The virtual identity of the two theories, notwithstanding 
their completely different formulations, can be considered 
as a source of insight for the understanding of the physical 
problem. The fact that the results can be obtained with the 
much simpler creep tide model, before the introduction of 
the elastic tide, is also insightful. 

2. Two angles play a major role in tide theory, the geode¬ 
tic lag (eo)i which measures the asymmetry of the shape of 
the body, and the phase of the semidiurnal argument in the 
anelastic tide ((Tq), which measures the asymmetry of the 
shape of the anelastic deformation. In rotating bodies not 
trapped in a spin-orbit resonance (like the Earth), with low 
relaxation factor, the geodetic lag is close to zero and the 
phase of the semidiurnal argument in the anelastic tide is 
close to 90 degrees. This result is the same in both theo¬ 
ries: the creep tide theory of Ferraz-Mello (2013) and the 
Maxwell model of Correia et al. (2014). It does not depend 
on the used theory. It is also true for modern versions of 
Darwin’s theory (Efroimsky, 2012a). 

Bodies trapped in a spin-orbit resonance were not con¬ 
sidered in this paper since, in that case, n ^ 0 and then the 
phase of the anelastic tide (ctq) tends to zero. This is also 
so in the frame of Efroimsky’s theories (Efroimsky, 2012b; 
Williams and Efroimsky, 2012). In the particular resonant 
case of Mercury, the equations are somewhat different be¬ 
cause the dominant tide has a different argument {v — n 
instead of v, see Ferraz-Mello 2015). The anelastic tide is 
then dominated by the term whose phase is 

V — n 

(T_i = arctan- > 0. 
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3. The Earth’s lag determined by Ray et al. (2001) from 
the combination of tidal solutions from satellite tracking 
data and from Topex/Poseidon satellite altimeter data is 
the geodetic lag. The same data may be used to determine 
the phase of the semidiurnal argument in the anelastic tide 
and gives cto = 89.80 ± 0.05 degrees. 

4. The results of the Maxwell model of Correia et al. 
(2014) can be obtained in a much simpler way with the 
methodology used in Ferraz-Mello (2013). The main results 
can be obtained from the solution of the Maxwell model’s 
creep (Eq. [T71) and a few elementary physics laws. The elas¬ 
tic component of the tide and its effects can be added af¬ 
terwards. 
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